Josephson Effect through an Isotropic Magnetic Molecule 
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We investigate the Josephson effect through a molecular quantum dot magnet connected to super- 
conducting leads. The molecule contains a magnetic atom, whose spin is assumed to be isotropic. It 
is coupled to the electron spin on the dot via exchange coupling. Using the numerical renormaliza- 
tion group method we calculate the Andreev levels and the supercurrent and examine intertwined 
effect of the exchange coupling, Kondo correlation, and superconductivity on the current. Exchange 
coupling typically suppresses the Kondo correlation so that the system undergoes a phase transition 
from to TT state as the modulus of exchange coupling increases. Antiferromagnetic coupling is 
found to drive exotic transitions: the reentrance to the tt state for a small superconducting gap 
and the restoration of state for large antiferromagnetic exchange coupling. We suggest that the 
asymmetric dependence of supercurrent on the exchange coupling could be used as to detect its sign 
in experiments. 

PACS numbers: 73.63.-b, 74.50.+r, 72.15.Qm, 73.63.Kv 
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Molecular spintronics [ij aims at exploring spin- 
dependent electronic transport through molecules with 
intrinsic degrees of freedom such as spin, connected to 
leads of various nature. On the theoretical and experi- 
mental side, recent advances have concerned both coher- 
ent and incoherent [1, 0, S] transport through these 
molecular quantum dot magnets (MQDM). They consist 
of a magnetic molecule with either a large or a small 
anisotropy, as is the case for a endofuUerene molecule Q ■ 

Here, we provide a nonpcrturbative computation of the 
low temperature transport properties of a MQDM con- 
nected to superconducting leads using a numeral renor- 
malization group (NRG) approach. The Josephson cur- 
rent allows a diagnosis of the interaction between the 
intrinsic spin of the molecule, its itinerant electron spin, 
and the polarization of the leads. It has been known for 
some time (1, 0, [lol 
periments |12| . 
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and recently analyzed in ex- 
that a quantum dot sandwiched between 
superconducting leads can show a tt junction behavior 
[l3|. At the same time, a quantum dot connected to 
leads at low enough temperatures exhibits the Kondo 
effect 14 1 . It was shown Iq, 15, 1^ that with supercon- 



ducting leads, at low temperature the junction state of 
the Josephson current is restored when the Kondo tem- 
perature exceeds the superconducting gap. The stability 
of this Kondo phase is put in question in the presence of 
additional spin degrees for freedom [l3| which may com- 
pete with Kondo screening. Here the Josephson current 
flows through an isotropic MQDM which can describe 
a endofuUerene molecule [l8[. The electron spin in the 
quantum dot and the magnetic ion inside it interact via 
an exchange coupling [4| . We calculate the Andreev level 
(AL) spectrum and the supercurrent and determine the 
spin of the ground state. We find that the exchange cou- 
pling typically suppresses the Kondo effect and drives a 
transition from to tt state. Moreover, antiferromag- 
netic coupling is found to drive exotic transitions: the 



reentrance to tt state for small superconducting gap and 
the restoration of state for large J . 

The MQDM connected to two s-wave superconducting 
leads (inset of Fig. [1]) is modeled by a single-impurity 
Anderson model: TL = TLm + Hl + Ti-r, where 

(1) 
(2) 



T^M = f'On + Un^rii + JS • Sc 



(3) 



Here c^k^i (d^) destroys an electron with energy ek, and 
spin fi on lead £ = L, R (on the carbon cell); n^k and n 
are occupation operators for the leads and the cell. The 
single-particle energy eo can be tuned by gate voltages. 
J denotes the exchange energy between the ion spin S 
and the electron spin So = '^Ji^'mm'^m'- ^ 

superconducting gap. Except for the finite phase dif- 
ference (j) = (j^L — (t>R, the leads are identical and their 
coupling to the MQDM is symmetric. The hybridiza- 
tion between the molecule and the leads is well char- 
acterized by a tunneling rate T — Trpolip, where po is 
the density of states of the leads at the Fermi energy. 
As we are interested in the low temperature behavior, 
we concentrate for the most part on the Kondo regime 
with a localized level — eo 3> F with large charging en- 
ergy U :s> \eo\. Specifically, we choose eo = —O.ID (the 
band width D is taken as a unit of energy), F = O.OII?, 
and U = oo and introduce the bare Kondo temperature 
T° = /Dl72exp (1 + f)] (at J = A = ^YThe 
energy spectrum is found with the NRG method [19] ex- 
tended to superconducting leads [l^, Within the 
NRG method, the supercurrent is directly obtained by 
evaluating the expectation value of the current operator 

Fig. [H shows the phase diagram of our system, which 
constitutes the main result. The jmiction property 



FIG. 1: (color online) Schematic phase diagram of a MQDM 
superconducting junction system [see the upper inset] indi- 
cating the 0, 0' (blue), tt' (green), and tt regions. Each region 
is divided into two subregions according to the ground-state 
S and S—1/2 for Oj'' and O2 regions and 5—1/2 and 
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5+1/2 for TT^'' and TTj regions, respectively. Note that the Oi 
state exists only along the line J=0 [see the lower inset]. For 
larger molecular spin S'>S (see the dotted lines), the phase 
boundaries between Oi and 7ri/2 are shifted toward smaller 
jjj, and one between O2 and tti moves toward larger J. 



switches between and tt state, depending on the 
strengths of J and A with respect to T^. For J = 0, 
the system undergoes the Kondo-driven phase transition 
B EBl '■ The ground-state wave function is of spin sin- 
glet kind for A < Ac « 1.84 T]^ and of spin doublet for 
A > Ac. In the strong coupling limit (A < Ac) Kondo 
correlations screen out the localized spin and Cooper 
pairs tunnel through the Kondo resonance state, result- 
ing in a 0-j unction In the weak coupling limit 
(A > Ac), strong superconductivity in the leads leaves 
the local spin unscreened and the tunneling of Cooper 
pairs subject to strong Coulomb interaction acquires an 
additional phase tt, making a 7r-junction 15, T^ . 
It is also found [15|| that the transition is ^-dependent so 
that a narrow region of the intermediate states 0' and tt' 
exists; see the enlarged view in Fig. [TJ 

Finite exchange coupling between electron spins and 
the ion spin introduces another electronic correlation and 
affects Cooper pair transport. Fig. [2] shows typical varia- 
tions of ALs and supercurrents with J along the line aa' 
(see Fig. [U in the strong couphng limit (A/T]^ = 0.1). 
Any finite J clearly induces a splitting in subgap exci- 
tations and consequently causes a crossing between the 
ground state and the lowest excitation at </) / tt (at least 
for jJ/T^I ^ 0(1)); the level crossing otherwise takes 
place only at (f> — tt. Across the crossing, the ground 
state spin is changed from S to 5*+ 1/2 for J ^ 0. Sim- 
ilarly, the ALs defined as the one- electron/ hole subgap 
excitations (identified as the poles of the dot Green's 
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FIG. 2: (color online) (a) ALs in units of A and (b) supercur- 
rents I in units of 1^ = eA/h as functions of <f) in the strong 
coupling limit {A/T^ =0.1) for various values of J/T^: see 
the hue aa' in Fig. [T] Here the ion spin 5* is set to 1/2. 



functions [2l|) exhibit discontinuities like kinks in the 
spectra; for J ^ two outmost ALs with spin S'±l/2 
with respect to the spin-S* ground state cannot remain 
as (spin-1/2) one-electron excitations with respect to the 
ground state with spin 5+1/2 at the transition and are 
replaced by new ALs with spin 5+1. In parallel with 
an abrupt change in ALs, the supcrcurrent-phase rela- 
tion (SPR) shows a discontinuous sign change (note that 
/ oc ~dEA/d6,as the continuum-excitation contribution 
is negligible [2l|), culminating in a transition from to tt 
state: two tt*^'^ states labeled as tt^^\ are identified accord- 
ing to the ground-state spin 5 + 1/2, respectively. The 
intermediate states 0'^ and tt'-^^^ defined as in Ref. 
111 ]. The full state exists only at J = because any 
small J drives the system to the tt state at = tt; see 
Fig. [2] The curve of /(</)) then has three distinct segments 
jl5j . The central segment resembles that of a short bal- 
listic junction, while the two surrounding segments are 
parts of TT-junction curve. As J grows in magnitude the 
central segment shrinks and eventually vanishes. The 
SPR then becomes sinusoidal like in a tunnel junction. 
It should be noted that the O-tt transition is asymmetric 
with respect to the sign of J: the transition for J > 
takes place at SEs ~ T^, where SEs = -5^(25+1) is the 
exchange-coupling energy gap, while the state survives 
much larger ferromagnetic coupling {J < 0). Once the tt- 
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junction is fully established, stronger ferromagnetic cou- 
pling does not lead to any qualitative change in the SPR, 
while a second transition back to state is observed for 
large antiferromagnetic coupling {J ^ A). The NRG re- 
sults distinguish the second state (O2) from the former 
one (Oi) in three points: (1) the ground state has spin 
5 — 1/2 like the tti phase, (2) the SPR is that of a tunnel- 
ing junction, and (3) the 7ri-02 transition has no interme- 
diate state. Figure [3] (c) shows that the critical current 
has its maximum at J = and decreases with increasing 
1 J| rapidly across the phase boundary for J > or rather 
gradually for J < 0. The critical current totally vanishes 
at the 7ri-02 boundary and increases again slowly with J 
in O2 phase (see the curve for A/T]^ = 0.01). 

The O-TT transitions (Oi-tti and 0i-7r2) can be attributed 
to the competition between superconducting and Kondo 
correlations as in the absence of exchange coupling. The 
relevant parameters are then the Kondo temperature Tk 
and the superconducting gap A, and the O-tt phase tran- 
sition occurs when they are comparable to each other: 
In our choice of parameters the transition happens at 
Ac/Tr- « 1.84. The exchange coupling manifests itself 
by renormalizing the Kondo temperature Tk{J)- To see 
this, we applied the poor man's scaling theory to a corre- 
sponding Kondo Hamiltonian with no superconductivity 
and 5 = 1/2: TLkm = J2k ^^kfik -|- JS • Sc 4- ( JkSb -I- JmS) • 
Sl, where Sl is the spin operator for the lead electrons 
at molecule site. The last term S • Sl describing direct 
coupling between spins of the ion and the lead electrons 
arises during the scaling process. The renormalization 
group analysis leads to the following scaling equations: 
together with J w J(A = D), 
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(4) 



As the band width A is decreased from D to Tk, the 
coefficient Jk, responsible for the Kondo correlation, di- 
verges and the scaling breaks down. In the presence of 
finite exchange coupling, however, since JJkJm > with 
Jm(A = I?) = and |Jm| ^ ^k, the term proportional 
to J in Eq. ([4]) turns out to slow down the flow of Jk and 
accordingly lowers the Kondo temperature. This point is 
confirmed by NRG calculations applied in the absence of 
superconductivity. As can be seen in Fig. [3] (a) and (b), 
the width of the spectral density for dot electrons, iden- 
tified as the Kondo temperature Tk{J), decreases with 
increasing | J| (for J < this decrease, being marginal, 
is not clearly shown with the logarithmic scale). We find 
out that for the ferromagnetic case the ratio Tk{J)/T^ 
coincides with Ac(J)/Ac(J = 0). For the antiferromag- 
netic case, the Kondo correlation is observed to be sup- 
pressed not only by the Kondo peak narrowing but by 
lowering the peak height. 

Antiferromagnetic exchange coupling can, on the other 
hand, exert a more profound effect than simply renormal- 
izing the Kondo temperature: it gives rise to a reentrant 
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FIG. 3: (color online) Spectral weights A{uj) for dot elec- 
trons coupled to normal leads with antiferromagnetic [(a)] 
and ferromagnetic [(b)] exchange coupling to ion spin for var- 
ious values of J/T^ (as annotated), (c) Critical currents as 
functions of J/T^ for different values of A/T]^ (see the anno- 
tations). The arrows locate transition points corresponding 
to data with the same color. Here we have used S = 1/2. 



transition to the tt state at small A and restoration of 
the state for large J. It is known that small anti- 
ferromagnetic exchange coupling (J < T]^), studied in 
the context of coupled impurities [22| and side-coupled 
quantum dot systems [2^ and observed in experiments 
24| . can produce a two-stage Kondo effect. After the 
magnetic moment of the dot is screened by conduction 
electrons below Tk, at a much lower energy scale (de- 
noted as T^) the ion spin is screened by the local Fermi 
liquid that is formed on the dot. T^ is then the Kondo 
temperature of a magnetic moment screened by electrons 
of a bandwidth ~ Tk and density of states ^ 1/{ttTk) 
[H: T^^ ~ rKexp[-^]. The second Kondo effect 



leads to a Fano resonance and makes a dip in the dot 
electron density of states as shown in Fig. [3] (a). The dip 
becomes widened with J and overrides the Kondo peak 
when Tk ~ Tk so that the Kondo effect is completely 
overridden. As long as A > Tj^, the second Kondo effect 
does not appear since the superconducting gap blocks 
any quasi-particle excitation with energy less than A. 
For A < TjJ, however. Cooper pairs notice the suppres- 
sion of the Kondo resonance level, and their tunneling 
is governed by cotunneling under strong Coulomb inter- 
action, forming a 7r-junction again. Since decreases 
with decreasing J, Ac decreases to zero as J ^ 0. Note 
that the extremely small TjJ < Tk (unless 6Es ^ T^) 
might make it hard to detect the reentrance even under 
rather weak thermal fluctuations with Tk > T > T^. 

The revival of the 0-state for strong antiferromagnetic 
coupling can be exp lained in the picture of cotunneling 
of Cooper pairs [lOj. In weak coupling limit, the fourth- 
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FIG. 4: (color online) (a) ALs in units of A and (b) super- 
currents I in units of as functions of 4> with J IT% = 10 
and IS./T% = 0.02 (at = —0.1) while the gate voltage td is 
tuned from -0.1 to -0.04. See the line hh' in Fig. [T] 



difficult to achieve, the relative strength J /T^^ can be 
controlled by the gate voltage which can tune the Kondo 
temperature. Fig. [3] proposes a possibility to observe a 
double transition (along the line hh' in Fig. [T]) as the gate 
voltage is swept. Note that the double transition is an 
evidence of strong exchange coupling (J » » A): 
for examples, with T]^ ~ 3K measured in a recent Cgo 
single- molecular transistor 2J|, one estimates J ~ 30K. 



Asymmetry of the phase diagram enables the sign and 
possibly the amplitude of J to be determined without 
ambiguity by observing the evolution of the SPR or the 
critical current. 

The authors thank W. Wernsdorfer, F. Balestro, and 
Mahn-Soo Choi for helpful discussions. This work is sup- 
ported by ANR-PNANO Contract MolSpintronics No. 
ANR-06-NANO-27. 



order perturbation theory leads to the supercurrent: 
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25-1-2 



2S+1 \Ek+Ek' SEs+Ek+Ei,, 



w here Ej, = ^A ^ + ej, Uk = V(T+ek/E^, Wk = 
V(l - ek/-Bk)/2, and fk = -e<i - jiS + 1) - E^ < 0. 
For antiferromagnetic coupling, the ground state for the 
uncoupled system has spin S — 1/2. After one electron 
tunnels through the molecule the system can be in spin 
eigenstate of either 5—1/2 and S + 1/2. The latter vir- 
tual process, costing more energy by the gap SEs, turns 
out to acquire a tt phase, contributing to a negative su- 
percurrent. The larger amplitude of this process by a 
factor 25-1-2 (degeneracy of the spin state 5-1-1/2) dom- 
inates over spin-preserving process as long as the gap 
6Es is small. For a large gap SEg, however, this process 
becomes negligible and the sign of the supercurrent is re- 
versed. Note that according to Eq. ^ the SPR is always 
sinusoidal and the current should vanish at the transi- 
tion, which is also confirmed in our NRG calculations. 

The physical arguments for the O-tt transitions dis- 
cussed so far are valid for arbitrary values of the ion spin 
5, while the phase boundaries are shifted with chang- 
ing 5 as shown in Fig. [1] The exchange-coupling energy 
gap SEg that is supposed to compete with Tk increases 
with 5 so that for larger 5 the transitions can occur at 
smaller J. On the other hand, we have observed that the 
7ri-02 transition takes place at slightly larger J for larger 
5. This is because the increase in the degeneracy factor 
25 -|- 2 overwhelms the decrease in matrix elements due 
to a larger energy cost by dEs [see Eq. 

Finally, we present potential experimental manifesta- 
tions of exchange-coupling-driven O-tt transition. While 
the direct control of exchange coupling in molecules is 
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